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ABSTRACT 



We have reanalyzed a data set of 99 low redshift (z < 0.1) Abell clusters 
studied previously by Rhee, van Haarlem & Katgert (1989), and determined 
their shapes. For this, three different measures are used: two of which were 
originally used by Rhee et al., and one of which was used by Plionis, Barrow 
& Frenk (1991) in their investigation of clusters in the Lick catalogue. We use 
Monte-Carlo simulations of clusters to investigate the errors in the methods. 
For low ellipticity, all methods overestimate the cluster elongation, whereas the 
opposite is true for a highly flattened system. Also background galaxies and 
shot noise have a rather large influence on the measured quantities. 
The corrected distribution of cluster ellipticities shows a peak at e ~ 0.4 and 
extends to e ~ 0.8, consistent with results of some previous studies. However, 
the present study uses more than twice the number of clusters as the earlier 
studies, and is self-consistent. That is, with the corrected distribution over 
projected cluster shapes we can reconstruct the observed distribution over pro- 
jected cluster shapes and the observed relation between the number of galaxies 
in a cluster and its ellipticity. To achieve this, we have to assume that there is 
an anti-correlation between the true (projected) ellipticity of a cluster and its 
number of galaxies. It is not necessary to assume that the ellipticity of a cluster 
increases when one only includes the brighter galaxies (as suggested by Binney 
1977). 

Using a redshift-independent richness criterion of Vink & Katgert (1994), it is 
shown that the richer clusters are intrinsically more nearly spherical than the 
poorer ones. Furthermore, the corrected distribution of cluster shapes is found 
to be more consistent with a population that consists of purely prolate clusters 
than with a purely oblate population. 

We compare the corrected true distribution of (projected) ellipticities with pre- 
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dictions from N-body simulations. For this, we use a catalogue of 75 N-body 
simulated clusters (van Kampen 1994) which assume a CDM spectrum with 
SI = 1.0. The simulations include a recipe for galaxy formation and merging. 
The model clusters are expected to be a representative sample of all real clus- 
ters. They show a good resemblance with the data both in radial profile and 
number of galaxies. 'Observing' these simulated clusters in exactly the same 
way as the real clusters produces an ellipticity distribution that extends to 
much higher e and that has too few nearly spherical clusters. Preliminary re- 
sults of simulations of the formation of clusters in an SI = 0.2 universe suggest 
that, on average, clusters are more nearly spherical in this case, as is expected 
on theoretical grounds. This shows that the elongations of clusters can provide 
a useful constraint on the value of il. 

Key words: Galaxies: clustering - cosmology: observations - large-scale struc- 
ture in the Universe. 

1 INTRODUCTION 

Since clusters of galaxies have only just collapsed or are still collapsing (e.g. Regos & Geller 
1989; van Haarlem 1992; Zabludoff & Franx 1993), their shapes will probably contain clues 
about the formation of large-scale structure and about the initial density field out of which 
structure in the Universe is believed to have formed by gravitational collapse. 

Most clusters, like elliptical galaxies, are non-spherical, and their shape is (as is the case 
for ellipticals) not due to rotation (Rood et al. 1972; Gregory & Tifft 1976; Dressier 1981). 
The perturbations that gave rise to the formation of galaxy clusters are likely to have been 
aspherical initially (Barrow & Silk 1981; Peacock & Heavens 1985; Bardeen et al. 1986). 
The asphericities are then amplified during gravitational collapse as the collapse proceeds 
more rapidly along the short axis (Lin, Mestel & Shu 1965; Icke 1973; Barrow & Silk 1981). 
In cases where there are good measurements, the elongation seems to be due to a velocity 
anisotropy of the galaxies (Aarseth & Binney 1978), as in elliptical galaxies. Binney & Silk 
(1979) claim that the elongation of clusters originates in the tidal distortion by neighboring 
protoclusters. Recently, Salvador-Sole & Solanes (1993) also suggested this by developing a 
simple model, with which they can explain the elongations of clusters. They point out that 
the main tidal distortion on a cluster is by the nearest neighboring cluster having more than 
45 galaxies. Their model also explains the (weak) alignment between neighboring clusters 
(Oort 1983; Binggeli 1982; Rhee & Katgert 1987; Plionis 1993) and the strong alignment 
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between clusters and their first ranked galaxy (Sastry 1968; Carter & Metcalfe 1980; Dressier 
1981; Binggeli 1982; Rliee & Katgert 1987; Tucker & Peterson 1988; van Kampen & Rliee 
1990; Lambas et al. 1990; West 1994). Recently, Bertscliinger & Jain (1993) and van de 
Weygaert & Babul (1994) claimed that shear fields probably have a large infiuence on the 
shapes of clusters. This shear can be due to the clusters own asphericity or can be induced by 
the surrounding mass distribution. It can break up a single cluster into more distinct ones, or 
promote a smaller cluster to a larger one. Another possible reason for the large elongations 
of clusters is given by Zeldovich (1978), who claims that it is due to the fact that clusters are 
born in sheets of gas ("pancakes"). However, this would require structures to have formed 
via a top-down or fragmentation scenario. In such a scenario clusters fragment out of larger 
structures and will have a certain elongation in the direction of the sheets. Although this 
scenario can also explain the alignment between neighboring clusters, it is probably ruled 
out (Peebles 1993) because there are clear examples of galaxies that are assembled before 
the systems in which they are now found. Also, the existence of high-redshift quasars is a 
problem for the top-down scenario. 

Some N-body simulations have shown that the shapes of dissipationless dark halos of 
clusters most likely will contain very little information about the fiuctuation spectrum giving 
rise to the perturbations (West, Dekel & Oemler 1989; Salvador-Sole 1993). On the other 
hand Quinn, Salmon & Zurek (1986), Efstathiou et al. (1988) and Crone, Evrard & Richstone 
(1994) demonstrated that there is a link between the final cluster structure and the initial 
spectrum of density fiuctuations. It is not clear how this is for the luminous part of the 
cluster, i.e. the galaxies (van Kampen 1994). 

A cosmological parameter that certainly infiuences the distribution of galaxies in a cluster 
is the density parameter 0. Evrard (1993) and Evrard et al. (1993) have shown, by comparing 
N-body simulations with X-ray images, that clusters in a low-0 universe are much more 
regular, spherically symmetric, and centrally condensed than clusters in an = 1 universe. 
This is because in a low-0 universe, structures will collapse at a larger redshift than in a 
high-0 universe (Maoz 1990) and thus have had more time to virialize and wipe out their 
asymmetries and substructures. Eurthermore, the timescale in a low-0 universe is (much) 
longer than in an Einstein-de Sitter universe (see e.g. Padmanabhan 1993, eq. [2.76]). Evrard 
(1993) and Evrard et al. (1993) conclude that the complex nature of real cluster X-ray 
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morphologies is better matched by an Einstein-de Sitter model than by = 0.2 models, 
even if they include a cosmological constant in the latter. 

However, the observational information on the distribution of cluster shapes is not unam- 
biguous. While Rhee, van Haarlem & Katgert (1989; hereafter RHK) find that most clusters 
are nearly spherical, with a peak at e ~ 0.15, Carter & Metcalfe (1980) , Binggeli (1982) 
and Plionis, Barrow & Frenk (1991; hereafter PBF) find them to be much more elongated, 
with most clusters having e ~ 0.5 (e denotes the projected ellipticity of a cluster, defined as 
one minus the projected axis ratio). 

In this paper we try to remove the inconsistency between the different determinations of 
the cluster shapes by comparing the methods used by the various authors. We concentrate 
on the samples used by RHK and PBF and apply the method of PBF to the data of RHK. 

The outline of the paper is as follows: in Section 2, the sample of clusters is described. In 
Section 3, we describe the methods used by RHK and PBF. Section 4 outlines the Monte- 
Carlo simulations, used to correct for the different systematic and random errors in the 
various methods. In Section 5 we describe the corrected results and investigate the relation 
between the elongation of clusters and their richness. In Section 6, we try to put constraints 
on the distribution over the intrinsic ellipticities of the clusters by assuming them to have 
a single triaxiality parameter. In Section 7 we use the N-body simulations of clusters, taken 
from a catalogue of model clusters by van Kampen (1994), with either b = 1.0,0 = 0.2 
OT b = 2.0,0 = 1.0. The shapes of the luminous part of the clusters in these models are 
compared with the data. In Section 8 we give a discussion and summarize the results. 

A Hubble constant Hq = lOOh km s~^ Mpc~^ is used throughout the paper. 

2 THE DATA SETS 

The cluster samples are fully described by the authors (RHK and PBF). For completeness, 
the main properties are reviewed here. 

2.1 RHK sample 

The RHK sample consists of 99 Abell clusters of Abell richness i? > 1, redshift z < 0.1, 
lO'' < a < 18^^, S < —25° and b > 30° (the original sample consists of 107 clusters, but only 
for the 99 clusters in Figure 5 of RHK all data were available). The sample is complete up to 
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a redshift z = 0.08. The galaxy positions and brightnesses were determined using the Leiden 
Astroscan automatic plate measuring machine, from glass copies of the red Palomar Sky 
Survey plates. The machine has been described by Swaans (1981). To identify the galaxies, 
the following stages were gone through: first, all objects down to a certain minimum threshold 
level in surface brightness were detected. Then stellar objects were separated from the non- 
stellar objects, using a separation criterion based on that devised by Kron (1980). The object 
position was then defined as the mean position of the five brightest pixels of the nine pixels 
surrounding the first object position. The galaxy catalogues have magnitude limits of niR 
between 18.0 and 19.0. 

The cluster centers were roughly determined using the coordinates in mm given by Sastry 
& Rood (1971). Then the brightest galaxy within ~ 0.25h~^ Mpc was found, and all galaxies 
which were detected within lh~^ Mpc from this brightest galaxy, were selected; note that 
RHK used h = 0.5 in their paper. Later, Vink & Katgert (1994) redetermined the centers 
by defining a circular area around the old center and calculating the center of mass of the 
galaxies in this area, which then became the new center. This procedure was repeated until 
the center did not change anymore. Convergence was always achieved. 

To estimate the background, 16 randomly chosen field regions were scanned in exactly 
the same way as the cluster regions. Based on the galaxy counts in these 16 regions, it was 
found that the average field contamination to mji ~ 19 is ~ 50% within lh~^ Mpc from the 
cluster center. 

2.2 PBF sample 

The sample used by PBF consists of ~ 100 clusters (with \b\ > 40°) in the Lick map, based 
on the Shane & Wirtanen galaxy catalogue (Shane & Wirtanen 1967; Seldner et al. 1977). 
This catalogue contains 810,000 galaxies with mj, < 18.8 binned in 10 X 10 arcmin^ cells. 
PBF smoothed the counts using a discrete convolution with the weights used by Shectman 
(1985). They then generated three cluster catalogues, using three overdensity thresholds viz. 
cr / < (T>= 3.6, 2.5 and 1.8, corresponding to 5.0, 3.5 and 2.5 galaxies per cell, respectively (cr 
denotes the number of galaxies per cell). They based their discussion of cluster ellipticities 
mainly on the 3.6 catalogue. In order to qualify as a cluster, a surface density peak needed to 
occupy at least 5 cells that are above the threshold (such cells were called 'active cells'). The 
clusters in the 3.6 catalogue have redshifts up to 2; ~ 0.15, with a median value of ~ 0.07. 
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3 THE METHODS 

In this section, we outline the methods used by RHK and PBF to determine the shape of a 
cluster. 

3.1 Description 

RHK used two different methods, one based on moments (from now on called the moments 
method, its result indicated by ei), the other based on the inertia tensor (henceforth called 
the tensor method, its result denoted by £2). 

For the moments method, the various second order moments of the galaxy distribution 
are defined as 

l^ab = ^Y.^^y' ' («,& = 0,I,2; a + b = 2) , (I) 

where the sum is over all galaxies, x and y are the coordinates of the galaxy relative to the 
cluster center, and Ngai is the number of galaxies in the cluster. The cluster ellipticity is 
then given by 



where 



e = ^ . (3) 

For the tensor method, the inertia tensor is defined as 

1^3 =Y.^'^3/^^ ■> (^^' = 1.2) , (4) 

where the sum is over all galaxies, Xi = x, X2 = y and r is the projected distance of a galaxy 
to the cluster center. The cluster ellipticity is then given by 

^2 1 ^min / -^max ? (^) 

where 



^V(/ll + /22)2- 4(722/22-/12/21) (6) 

Both methods used by RHK calculate e within I Mpc from the center. A few clusters 
are located near the edge of a plate and e is then evaluated at a smaller radius. 
For the PBF method, the inertia tensor is defined as 
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Figure 1. Original distributions over ellipticities of RHK and PBF. The solid line gives the distribution for the moments 
method of RHK, the dotted line denotes the tensor method of RHK (both derived from the RHK data), and the dashed line 
indicates the results from the tensor method of PBF and as derived from the PBF data. 



where the sum is over all active cells having more than a = 3.6 < a >= 5.0 galaxies, 
Xi = X and X2 = y are the coordinates of the center of the cell with respect to the cluster 
center, and m is the number of galaxies in that cell (after smoothing). The cluster ellipticity 
is then defined according to eqs. (5) and (6), and will be denoted by £3. Only clusters having 
at least five active cells are included. This causes a systematic selection against non-centrally 
concentrated clusters. 

3.2 Original results 

The apparent (i.e., as originally measured by the authors) distributions over ei, €2 and £3 are 
shown in Figure 1. The solid line corresponds to the moments method of RHK, the dotted 
line to the tensor method of RHK, and the dashed line to the tensor method of PBF. It is 
obvious from Figure 1 that the two methods of RHK give quite different results, especially 
at low €apparent- The tcusor mcthod shows a lack of spherical clusters relative to the moments 
method. Since they are applied to exactly the same data, this clearly illustrates the effect of 
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using different methods. The PBF method is even more different and finds many eiongated 
ciusters. However, this may partiy be due to the different iimiting magnitude of both sampies 
(see Section 2). The figure shows the apparent inconsistency between different methods and 
maybe different data sets. 

Because the RHK sampie is aimost fuiiy avaiiabie to us, from now on oniy this data set 
wiii be used and the presented resuits for the PBF method are obtained by appiying this 
method to the RHK data set rather than to the PBF data set. Since the Lick; map oniy 
contains objects with a biue apparent magnitude nib < 18.8 and the RHK data set is based 
on red piates, one has to talce a corresponding iimit in the red band to compare both data 
sets. We have determined this iimit by ioolcing at the cumuiative number of active ceiis (of 
the RHK sampie) as a function of limiting magnitude (in the red), and compared this to 
the distribution over active cells of PBF (their Figure 7). The magnitude limit that gives 
the best fit is niR ~ f7.0. This means, that if we apply a magnitude limit niR ~ f7.0 to 
the RHK sample, we on average select the same galaxies as does the Lick catalogue. The 
conversion from blue to red magnitudes is in agreement with our photometric calibration 
(Vink & Katgert f 994) and with the average color B — R f .8, given by Carter & Metcalfe 
(f980) and Postman, Geller & Huchra (f988). Furthermore, we assume the cluster centers of 
the RHK sample to be randomly positioned with respect to the grid of cells that is needed 
to apply the PBF method. This may introduce some errors, because the measured ellipticity 
of a cluster is found to depend on the positioning of the center with respect to the cells. As 
mentioned earlier, a rather large fraction of the clusters are rejected because they do not 
have the 5 active cells required to include them in the calculations. 

4 MONTE-CARLO SIMULATIONS 

In this section we apply the methods described in the Section 3.1 to computer-generated 
(Monte-Carlo) clusters. In this way the relations between input ('true') and output ('appar- 
ent') ellipticity can be established. These relations are then used to statistically correct the 
distributions of measured ellipticities for systematic effects in the methods. 

4.1 Generating synthetic clusters 

For the model clusters we assume a triaxial density distribution according to a modified 
Hubble-profile, and stratified on similar, concentric ellipsoids. 
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p{m) = po[l + m^]-^/^ , (8) 
where 

x'^ 

'n' = -; + ^ + -, , a>b>c , (9) 

(Sarazin 1986, and references therein; RHK). po is the central galaxy density in the cluster. 
This profile, together with assumed values b/a and c/a of the intrinsic ellipticities, serves as 
the probability function from which the galaxies are drawn randomly. The projected surface 
density profile, associated with equation (8), is 

pirn) = po[l + m^]~^ , (10) 
where po = 2a po and 



™' = ^ + |j. (11) 

with x^y the projected coordinates on the sky, and a the cluster core radius. A value a ~ 
0.25h~^ Mpc is adopted for the core radius (Sastry & Rood 1971; Bahcall 1975; Dressier 
1978; Binggeh 1982; Jones & Forman 1984; Sarazin 1986, and references therein; RHK). See, 
however. Beers & Tonry (1986), who claim that the core radius may be much smaller or 
even disappear if not the galaxy positions but e.g. cD galaxies are used to determine the 
cluster center. 

Furthermore, we include a background in the simulations. The background galaxies are 
distributed at random over an area of [—2,2] X [— 2,2]/i~^ Mpc around the cluster center. 
Such a large area is chosen because the PBF-method is not limited to a fixed area, so one 
has to include enough background. To include the right number of background galaxies, the 
magnitude distribution of the field galaxies in the 16 sky fields of RHK is used together 
with an assumed magnitude limit. These field galaxy counts are in agreement with results 
of Binggeli (1982), Jones et al. (1991) and Metcalfe et al. (1991). Although the background 
galaxies will, in reality, not be distributed randomly over the area because galaxies are 
believed to lie in walls and filaments, and those that are near a cluster may experience the 
infiuence from that cluster, the random positioning is found to give a good approximation. 
In fact, the results depend only very slightly on the precise number of background galaxies. 

Similar simulations are done for core radii of 0.1 and 0.5h~^ Mpc. They give very similar 
results for the correction cubes, and so our results will be rather insensitive to the actual 
value of the adopted core radius. 
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Figure 2. Example of a Monte-Carlo cluster. This cluster has a true projected ellipticity of 0.02, contains 150 member galaxies, 
and 60 background galaxies randomly distributed over the area of [-2, 2] X [-2, 2]/i-l Mpc around the cluster center. The upper 
lefthand plot shows the galaxy distribution of one of the realizations. Circles are cluster galaxies, crosses indicate background 
galaxies. The upper righthand plot shows the variation of ei with radius, as derived from 500 realizations. Circles and solid lines 
are cumulative values, whereas crosses and dotted lines are differential values. The middle lines indicate the average value of the 
ellipticity over the 500 realizations. The upper and lower lines denote the maximal and minimal value of ei , respectively. The 
lower left plot shows the variation of £2 with radius. The lower right plot shows the variation of £3 with overdensity threshold. 

4.2 Runs 

The Monte-Carlo simulations are performed as follows: We construct clusters with two in- 
trinsic axial ratios b/a and c/a and for given viewing angles one can then calculate the 
projected ellipticity b/a. For a fixed projected ellipticity etme = b/a the galaxy distribution 
on the sky will always be the same, independent of the intrinsic three-dimensional shape. 
This is because of the assumption of similar, concentric ellipsoids (e.g. Stark 1977). For each 
projected ellipticity, etme = 0.00, 0.05, 0.10, 1.00, 500 clusters are generated, according 
to the above constraints. An example of such a Monte-Carlo cluster is shown in Figure 2, 
together with the ellipticities as measured by the three methods. This cluster has a true 
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projected ellipticity of 0.02, contains 150 member galaxies distributed according to equation 
(10), and 60 background galaxies randomly distributed over the area of [—2,2] X [— 2,2]/i~^ 
Mpc around the cluster center. The upper lefthand plot shows the galaxy distribution of one 
of the realizations. Circles are cluster galaxies, crosses indicate background galaxies. The 
upper righthand plot shows the variation of ei with radius. Circles and solid lines are cumu- 
lative values, crosses and dotted lines are differential values. The middle lines indicate the 
average value of the ellipticity over the 500 realizations. The upper and lower lines denote the 
maximal and minimal value of ei, respectively. The lower lefthand plot shows the variation 
of €2 with radius. The symbols and lines have the same meaning as in the upper righthand 
plot. The lower righthand plot shows the variation of £3 with threshold. The middle line 
again denotes the mean, the upper and lower lines denote the maximal and minimal values 
of £3, respectively. This figure shows some striking features. First, the apparent ellipticity 
changes with radius. This is best seen in the lower left plot. It is only caused by noise in 
the method because the cluster is constructed to have a constant etme with radius. Second, 
there is a large spread in apparent ellipticities, which decreases with radius. This is due to 
the smaller number of galaxies in the center. Third, the average measured tapparent is always 
larger than the true one (in this particular case of etme = 0.02; but see later). 

The above procedure is performed for clusters containing different numbers of member 
galaxies, Ngai, ranging from 50 to 1000. This generates a 'correction cube', which gives the 
distribution over eapparent for a range of etrue-values and a range of A^^a^-values. With this 
cube we are able to correct statistically for the systematic errors and noise in the different 
methods. Figure 3 shows a slice through this cube, for clusters having 100 member galaxies 
and 60 background galaxies. The contours are equidistant and refiect the probability that a 
true ellipticity etme results in a measured ellipticity [i = 1,2,3). 

4.3 Results of simulations 

For small numbers of galaxies , the moments method of RHK gives ei ~ 0.10 for etme ^ 0.15. 
For higher etme^ ^1 increases too, but is always well below the true value. At etme = 1-0, 
ei ~ 0.7. The method performs better for larger numbers of galaxies, i.e., the measured 
curve approaches the theoretical curve and the distribution over ei gets narrower, indicating 
that the spread in the measured values decreases. 

The results for the tensor method of RHK are similar to those for the moments method, 
except that the method is slightly less accurate for small etme (in absolute sense; the width 
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Figure 3. Relation between true and apparent ellipticity for clusters containing 100 member galaxies. A value of 60 background 
galaxies is adopted. The upper plot is for the moments method of RHK, the middle plot is for the tensor method of RHK and 
the lower plot is for the tensor method of PBF. The contours are at equidistant values. 
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Figure 4. Resulting distributions over etrue, as calculated by the Lucy iteration, and only showing the signal for '^true,mtn < 
etrue < '^true,max • A limiting magnitude of 17.0 is adopted. The lefthand figure is for the moments method of RHK, the 
righthand figure for the tensor method of RHK. 

of the distributions are similar), and more accurate for higher etme- The dependence on the 
number of galaxies is similar. 

The tensor method of PBF behaves differently from the other two. Especially for small 
numbers of galaxies per cluster and small etme^ the method is noisy. This is clear from Figure 
3, where almost no values of £3 between 0.1 and 0.2 are measured due to discretization effects. 
For larger etme and Ngai the method performs better, and for high Ngai it is similar to the 
tensor method of RHK. This is to be expected, since both methods basically are the same, 
except that in the method of PBF the galaxies are binned in cells. For large numbers of 
galaxies, the discretization due to the binning of the galaxies in cells becomes less significant 
since more cells are taken into account. 

The systematic biases in the three methods can be understood as follows. First, for 
small etrue, the average tapparent will always be too large since the problem is asymmetric. 
For example, for etme = 0.0, all values of eapparent will be larger than or equal to 0.0 and thus 
the average will also be larger than or equal to 0.0. This effect is analogous to that proposed 
by Thuan & Gott (1977) to account for the difference between the observed and predicted 
distributions of axial ratios in elliptical galaxies. For larger etme this problem becomes less 
significant. Obviously, for etme = 1-0, a similar problem arises, but now one expects too low 
an average value. Second, the first two methods use a circular aperture, which tends to make 
the clusters more spherical. This effect is more important for the more elongated clusters 
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and explains why the measured values are always too low for high etme- For the third method 
this does not hold, but this method is expected to suffer a lot from discretization noise, since 
most of the clusters have only 6-8 active cells (PBF; Plionis 1993). This explains the rather 
discrete behavior of £3 (Figure 3). Third, the background galaxies also tend to lower the 
values of tapparent^ siucc they are assumed to be distributed at random over the aperture. 
The effect will be largest for higher etme and will have less effect on £2 than on ei and £3, 
since in the former method the weight of a galaxy decreases with radius (see Section 3.1). 

The PBF method shows some additional features. First, a rather large fraction of the 
simulated clusters is rejected because they have less than five 'active' cells with five or more 
galaxies (corresponding to an overdensity of 3.6). Second, almost all rejected clusters have 
£3 < 0.1. There may even be an effect depending on redshift. At larger redshift, a cluster 
will in general occupy fewer cells. Therefore, it has a larger chance of being rejected. One 
may conclude that the spherical clusters are mainly rejected at larger redshift. This effect 
may be visible in Figure 15(a) of PBF, where the authors find a correlation between redshift 
and ellipticity, in that the more distant clusters have a larger elongation. They almost find 
no spherical clusters at larger redshifts. 

Note that Fasano et al. (1993) have applied a very similar procedure to determine the 
shapes of galaxy groups. However, they did not include background galaxies in the sim- 
ulations that investigate the main errors and biases. Furthermore, the average number of 
galaxies in groups is only about 10. This will increase the noise in the determination of the 
shape parameter considerably, as pointed out before. 

In the following section the correction cubes will be applied to the distributions of ap- 
parent ellipticities (Figure 1). 



5 CORRECTING THE DISTRIBUTIONS 

5.1 Lucy's method 

What we have is an observed distribution N{eapparent) over projected ellipticities and a prob- 
ability function p{eapparent\^true) giviug the probability of 'observing' tapparent if the cluster's 
true ellipticity is etme- We now want to deduce the distribution of true ellipticities, N{etrue)^ 
from this. This distribution is given by 

{^apparent) {^true)p{^apparent\^true) • i,^^) 
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Figure 5. Cumulative number of clusters, having a true ellipticity larger than etrue- This profile is obtained using the histograms 
in Figure 4. The solid and the dashed line are for the moments method of RHK at a limiting magnitude of 17.0 and 18.0, 
respectively. The dotted and dot-dashed line are for the tensor method of RHK at a limiting magnitude of 17.0 and 18.0, 
respectively. Both the total number of clusters in the sample, Ntot, and the total number of clusters with etrue < '^true,mtn 
are indicated. 

However, there is one additional problem: p{€apparent\^true) depends on the number of galaxies 
in a cluster. So what one has to do is to bin N{eapparent) according to the number of galaxies 
per cluster, Ngai- Then, for each A^^a^-bin, one has a N{eapparent)^ and one has to solve 
eq. (12). The standard way to do this is using the Lucy (1974) algorithm. This method 
guarantees smooth, positive solutions. However, the inversion from N{eapparent) to N{etrue) 
will be non-unique, as is obvious from Figure 3, because dcapparentl detme ~ for etme smaller 
than some critical value €true,min- For etme < ^true,mim all N{etrue) give the same contribution 
to N{eapparent)^ SO ouc cau uot discriminate between different A^(etrMe)-pi'ofiles in this etme- 
range. The same effect operates near etme ~ 1-0 because there also deapparent/ detme can 
become zero. If this happens, there is a critical etrue,max7 which is in general close to 1.0. 
Also, only few clusters are expected to have such large etme- Note that the values of etrue,min 
and etrue,max depend on the magnitude limit one adopts, the method used, and the A^^a^-bin. 
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Figure 6. Distribution over true projected ellipticities. 

Table 1. Number of clusters that have either etrue < '^true,mtn or etrue > '^true,max- A limiting magnitude of 17.0 is adopted. 
The total number of clusters in the sample is 99, though method 3 rejects more than one-third or these because they occupy 
too few active cells. 

method 1 method 2 method 3 



0.4 2.64 
0.3 32.26 
0.2 0.11 



^true,min ^min 

0.4 10.60 
0.3 18.47 



^true,min ^min 

1.0 57.86 
0.2 0.99 



^true,max ^max ^true,max ^max ^true,max ^max 

0.9 2.47 

For each Ngai-hin we do the Lucy-iteration, and after convergence one is left with a 
distribution over etme in the range €true,min ^ ^tme ^ ^true,max7 with the total number of 
clusters having etme < €-true,nitn given by 



"^true.min 
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Figure 7. EUipticity versus the number of galaxies in a cluster. The ellipticity is the apparent ellipticity as measured by the 
tensor method of RHK and as applied to the original data. A limiting magnitude of 17.0 was adopted. 

and with the total number ol clusters having etme > ^true,max given by 

yl.O 



^tru 



e.max 



Table 1 gives Smin and S^ax for all three methods and adopting a limiting magnitude of 
17.0. Srain and Sraax generally increase when adopting a larger limiting magnitude. This is 
because the number of background galaxies roughly increases by 10°'^™'"" (Binggeli 1982), 
i.e., a factor of ~ 4 per magnitude, while the number of cluster galaxies increases less 
rapidly. Then, we combine all A^^a^-bins by simply summing all N{etrue) in the intervals 
^true,min ^ ^true ^ ^true,max- The rcsults of this Summation are shown in Figure 4. The 
lefthand panel is for the moments method of RHK, the righthand panel for the tensor 
method of RHK. These diagrams thus show the contribution to the distribution over true 
projected ellipticities, resulting from the Lucy-iteration in the intervals where etme is uniquely 
determined. For the PBF-method, no diagram is shown, because this would be completely 
empty. That is, at a limiting magnitude of 17.0, the number of galaxies in the clusters is so 
small that the PBF-method has no intervals in etme where it can be uniquely determined. 

Now one has to take into account the number of clusters Smin and Smax- For each A^^a^-bin, 
we have the total number of clusters that have etme < ^tme mim or that have etme > ^tme max- 
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What one can do is plot the number ol clusters that have a true ellipticity larger than a 
certain etme- This cumulative distribution is shown in Figure 5, for both methods of RHK, 
and at limiting magnitudes of 17.0 and 18.0. The solid and dashed lines show the results 
for the moments method at limiting magnitudes 17.0 and 18.0, respectively. The dotted 
and dot-dashed lines show the results for the tensor method at limiting magnitudes 17.0 
and 18.0, respectively. As Smin is only non-zero below a certain €true,min (which depends on 
limiting magnitude and method), all profiles should coincide above this value. For example, 
at niiim = 17.0, both methods only have a non-zero contribution from Smin below etme = 0.4 
(the number Smax = 2.47 for the first method has already been accounted for in Figure 
5). Therefore, for etme > 0.4 the solid and dotted line should coincide. This is not exactly 
the case, but it is very well possible to select a best-fitting profile to all curves. Note that 
the dashed line (moments method, ninm = 18.0) shows that the number of clusters in all 
(etrue.mm, etrMe,maa;)-intervals drastically decreases when the limiting magnitude is increased 
(as stated before). This curve hardly puts any constraints on the total distribution of clusters 
over €true- Then, for a specific method and a certain niura, all contributions of Smin have to be 
added to the profiles in Figure 5, under the constraints that the curves should be consistent 
and that the total number of clusters is 99 (as indicated in the figure). This 'extrapolation' 
is performed with a very low order polynomial fit. That is, we assume there to be no rapid 
variations in the number of clusters at low etme- That it is indeed possible to make all curves 
nearly consistent guarantees that one gets a nearly unique solution. 

5.2 Results 

The resulting distribution over etme is shown in Figure 6, and is just the derivative of the 
'extrapolated' profiles in Figure 5. It looks rather different from the apparent distributions 
of RHK in Figure 1. Not only does the peak shift to e ~ 0.4, but the profile also extends 
to e ~ 0.8, whereas in Figure 1 there were no clusters with elongations larger than e ~ 0.4. 
This suggests that clusters are more elongated than the RHK-methods in Figure 1 would 
indicate. The PBF method underestimates the ellipticities much less, but this is partly due 
to discretization. 

To check the results, we construct 99 clusters satisfying both the profile of Figure 6 and 
the distribution over number of galaxies per cluster as in the original RHK data set. Then 
we apply all three methods to these data and compare them to the original observations. To 
include a possible relation between the ellipticity and number of galaxies of a cluster, the 
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Table 2. Assumed anti- correlation between cluster elongation and its number of galaxies for the different limiting magnitudes. 

^true ^gal ^gal ^gal 

{mirm = 17.0) {mu,m = 18.0) {mu,m = 19.0) 

[0.0,0.2] [150,500] [200,1000] [200,1000] 

[0.2,0.4] [50,200] [100,300] [100,350] 

[0.4,1.0] [0,100] [0,150] [0,200] 

apparent ellipticities (determined with the tensor method of RHK and for niura = 17.0) are 
plotted versus the number of galaxies in the cluster (Figure 7). Although the uncorrected 
ellipticities in Figure 7 will not be correct quantitatively, the following can be concluded: 
because for clusters containing many galaxies, e is rather well determined, the fact that 
£2 < 0.3 for Ngai > 150 indicates that there are no clusters in this A^^a^-range which have 
ttrue larger than 0.45 (as obtained from the correction cube). The clusters with fewer than 
~ 150 galaxies show measured ellipticities up to ~ 0.8. Because the true ellipticity is almost 
always larger than the apparent one, and because there really are clusters with etme > 0.45 
(see Figure 6), this means that these clusters (with Ngai < 150) must have rather large true 
ellipticities. Thus, there is a clear anti-correlation between cluster ellipticity and its number 
of galaxies (see also Section 5.3). The same anti-correlation is also found using the apparent 
ellipticities of the other two methods and limiting magnitudes of 18.0 and 19.0. For the 
method of PBF, the anti-correlation weakens towards lower limiting magnitudes, and for a 
limiting magnitude of 17.0 it is almost completely absent, as noted by PBF in their Figure 
15b. But this is due to noise and discretization, as discussed in Section 4.3. 

The anti-correlation we used in assigning an ellipticity and number of galaxies to the 
Monte-Carlo clusters is listed in Table 2. The Monte-Carlo clusters are constructed according 
to the result in Figure 6 and the distribution over number of galaxies per cluster from the 
original RHK data. The apparent elongations can be compared with the original data. 

The results are shown in Figures 8 and 9. In Figure 8 the lefthand column shows the 
profiles as obtained directly from the data set of RHK. The righthand column shows the 
^apparent distributions as obtained from the simulations. The limiting magnitudes are 17.0, 
18.0, and 19.0 (from top to bottom). Again, the solid line is for the moments method of 
RHK, the dotted line for the tensor method of RHK, and the dashed line for the tensor 
method of PBF. Note that the lower lefthand plot is somewhat different from Figure 1 
because the centers of the clusters were redetermined (see Section 2.1). The agreement 
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Figure 8. Results of simulating the data. The lefthand column shows the observed profiles over ellipticity, as calculated directly 
from the RHK data set (using the uncorrected ellipticities). The righthand column shows the 'observed' profiles over ellipticity, 
as reconstructed from the true profile of Figure 5 and from the distribution over numbers of galaxies. From top to bottom, the 
plots show the results for limiting magnitudes of 17.0, 18.0, and 19.0, respectively. The solid line is for the moments method of 
RHK, the dotted line is for the tensor method of RHK, and the dashed line is for the tensor method of PBF. 
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Figure 9. EUipticity versus the number of galaxies in a cluster. The ellipticity is the apparent ellipticity as measured by the 
tensor method of RHK and applied to the model clusters. A limiting magnitude of 17.0 was adopted. 

between real observations and the predicted 'observations' is good lor both RHK methods 
and at all three magnitude limits. For the PBF method the agreement is somewhat worse. 
The Kolmogorov-Smirnov probability for the observations and predictions to be the same 
is always larger than 0.97, except for £3 at a limiting magnitude of 17.0 (KS probability is 
0.83) and at a limiting magnitude of 19.0 (KS probability is 0.00). This, again, illustrates 
that the PBF method is unable to reproduce the original results because of the rather large 
scatter in the method. Figure 9 can be compared with Figure 7. The agreement is very good. 
Although the assumed anti-correlation described in Table 2 is rather simple, it is necessary 
to reproduce both the distributions over tapparent and the correlation between tapparent and 
Ngai in Figure 7. This indicates that the obtained profile over true projected ellipticities 
in Figure 6, together with the assumed anti-correlation between ellipticity and number of 
galaxies of Table 2, are at least consistent with the observations. Furthermore, it is clear that 
a etrue-distribution, which does not depend on magnitude, as is assumed in constructing the 
distribution over etme (see Section 5.1), can also explain the differences in observed profiles 
for different niura, i.e., while the observations show that the average ellipticity of a cluster 
increases for the more luminous galaxies, this may only be an artifact of the methods and due 
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to the smaller number of galaxies included. The upper row of Figure 8 also show a property of 
the PBF-method, mentioned previously. For lower limiting magnitudes this method rejects 
almost half of the clusters because they do not occupy 5 active cells. Also, it shows the 
discretization noise of this method: for high limiting magnitudes the result is similar to that 
of the tensor method of RHK, but for lower magnitudes it gets inconsistent and shows an 
almost flat distribution. 

The errors in Figure 8 are rather small and are obtained by trying different distributions 
over €true ^ud different anti-correlation schemes which all produce the same profile over 
^apparent ^s Originally observed, and which show the same relation between ellipticity and 
number of galaxies (see Figure 7). The errors are largest in the lowest bins. This is to be 
expected, since at low etme^ the measured quantity will be always very similar (see Figure 



Our results can be compared with earlier work of e.g. Carter & Metcalfe (f980) and 
Binggeli (f982). Carter & Metcalfe showed clusters to have an average ellipticity of 0.44 it 
0.04, somewhat larger than ours. However, their sample consisted of only 2f clusters and they 
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Figure 10. Relation between the richness of a cluster, as defined by Vink &z Katgert (1994), and its number of galaxies within 
lh~^ Mpc (evaluated at a limiting magnitude of 17.0). 
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Figure 11. Correlation between the apparent ellipticity of a cluster and its (redshift-independent) richness. The former is 
measured by the tensor method of RHK and for a limiting magnitude of 17.0. 

at which it is evaluated. The authors indeed found a decrease in eUipticity outwards. This 
was also noted by Flin (1984) and Burgett (1982). More inwards, there is a larger scatter 
due to the smaller number of galaxies involved. Binggeli (1982) showed the clusters to have 
high elongations with a peak at e ~ 0.45 and extending to 0.8. However, he only used 44 
clusters each with only 50 galaxies, thereby increasing the noise considerably. Thus, although 
the previous studies are consistent with our results, this study has much better statistics 
and shows internal consistency: we are able to reproduce all features in the observations: 
an apparent increase in ellipticity for the more luminous galaxies, and an anti-correlation 
between ellipticity and number of galaxies in the cluster. 

5.3 Relation between e and richness 

We have argued before (see Section 5.2) that the clusters containing more galaxies tend to 
be more spherical. However, an obvious bias is observed, in that the clusters containing the 
most galaxies, are most nearby. So it is not clear whether these clusters are intrinsically 
rich, or whether they just appear rich. To investigate this, we plot the richness of a cluster 
versus its number of galaxies within lh~^ Mpc and for a limiting magnitude of 17.0 (Figure 
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oblate spheroids (T = 0) 

Figure 12. Plane of intrinsic ellipticities {emt,l 5 '^tnt,2)- Only shapes with emt,2 ^ '^tnt,l are possible. Dotted lines are lines of 
constant triaxiality T. Oblate and prolate spheroids correspond to T = and T = 1, respectively. 

10). Here the richness is defined by Vink & Katgert (1994) as the number of cluster galaxies 
within the absolute magnitude interval [-22.5,-21.5]. This number is corrected for background 
contamination. This definition of richness makes it independent of redshift, in contrast to 
the original definition of Abell (1958). From Figure 10 it is obvious that there is a clear 
correlation between the richness of a cluster and its number of galaxies. Clusters of fixed 
redshift lie approximately on a straight line through the origin. For increasing redshift this 
line steepens in the direction of the arrow. For fixed richness, clusters at lower redshift have 
a higher number of galaxies. 

Together with the result of Figures 7 and 9 this means that rich clusters are intrinsically 
more spherical than poor ones (as we have assumed in simulating the observations in Section 
5.2). This is shown in Figure 11 where we plot the apparent ellipticity of a cluster versus its 
(redshift independent) richness. The correlation is even more clear than in Figure 7. Fven 
though the values of €2 will not be correct quantitatively, as in Figure 7, qualitatively the 
trend of this figure is correct. The same correlation was found by Struble & Ftaclas (1994) 
who calculated the shapes of 350 clusters of galaxies as a function of their shape and velocity 
dispersion. 

The physical reason for this correlation may be that regions of higher density turn around 
earlier from the Hubble fiow than lower density regions. Since we have a fixed aperture of 
lh~^ Mpc, richer clusters will have a higher overdensity inside this aperture. The mean 
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Figure 13. Inverted distribution of intrinsic ellipticities for three different values of the triaxiahty parameter T . The sohd hne 
is for T = 0.0 (oblate), the dotted line is for T = 0.5 (intermediate triaxial), and the dashed line is for T = 1.0 (prolate). 



relaxation time will also be shorter and this will cause the more spherical shape for rich 
clusters. This result is in disagreement with the result of PBF, who found no correlation 
between the shape of a cluster and its number of galaxies (their Figure 15a). 



Since the projected ellipticity results from the projection of two intrinsic ellipticities, the 
problem of determining the distribution over the two intrinsic ellipticities is undercon- 
strained. One possible way out is to assume all clusters to have the same triaxiality pa- 
rameter, r, defined as 



where a, 6, c are the axes of the clusters (see eq. [8]) and €int,i = 1 — (c/a) and €int,2 = 1 — (b/a) 
the intrinsic ellipticities, with €int,i the larger of the two. This procedure was also used by 
Franx, lUingworth & de Zeeuw (1991), who applied it to elliptical galaxies. 

Figure 12 shows the plane of all possible shapes {€int,ifynt,2)- Curves of constant T are 
shown. The limiting cases T = and T = 1 correspond to oblate and prolate shapes, respec- 
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Figure 14. Recovered distribution of true ellipticities for three different values of the triaxiahty parameter T . The sohd hne 
is for T = 0.0 (oblate), the dotted line is for T = 0.5 (intermediate triaxial), and the dashed line is for T = 1.0 (prolate). 

tively. We calculate the distribution over intrinsic ellipticities with Lucy's (1974) method, 
assuming all clusters to have the same triaxiality parameter T. This method guarantees 
smooth, positive solutions. The resulting distributions over €int,i are shown in Figure 13 
for T = 0.0 (oblate models, solid line), 0.5 (intermediate triaxial models, dotted line), and 
1.0 (prolate models, dashed line). All the distributions peak at €int,i ~ 0.55. In the oblate 
case, no clusters with €int,i ~ 0.0 occur, and a rather large percentage has a very large 
elongation. As expected, the purely prolate distribution extends to smaller €int,i and shows 
less very elongated clusters. This is because a prolate cluster has a rather small probabil- 
ity to be viewed as being spherical, so many of them are needed to produce the signal at 
small projected ellipticities. Oblate clusters have a rather small chance of being observed 
very elongated, so many of these are needed to produce the signal at very large projected 
ellipticities. The T = 0.5 case is intermediate. 

The Lucy iteration always gives a best fitting solution, and one does not yet know how 
good this fit is. To check this, we generate Monte-Carlo clusters according to the profiles in 
Figure 13 and then reconstruct the distribution over etme from this. The results are shown 
in Figure 14 and should be compared with Figure 6. For large elongations, all T-values are 



The Shapes of Galaxy Clusters 27 



50 ^ 
40 : 



^ 30 : 

o - 

o : 



10 



f"":"\ , I , , , I , I , frm , , , i 

0.0 0.2 0.4 0.6 0.8 1.0 

^apporent 

Figure 15. Distribution over observed ellipticities for clusters simulated in an Q = 1.0 universe. The lines have the same 
meaning as in Figures 1 and 7. 

equally good. However, the oblate distribution produces too many (nearly) spherical clusters. 
The other two distributions perform equally well. It is therefore to be expected that most 
of the clusters will have a triaxiality parameter between ~ 0.5 to 1.0. 

The following has to be remarked with regard to the intrinsic distributions. Both for a 
purely prolate and a purely oblate distribution there is a bias in the projected distribution 
because Abell clusters are defined as an overdensity in the surface density of galaxies, For 
prolate clusters, a cluster with its long axis pointing towards us has a larger probability of 
being detected since its surface density will in general be higher than that of a cluster with 
its short axis pointing along the line of sight. This means that for prolate clusters there 
is a bias towards observing more spherical than elongated ones. In this case, the projected 
distribution has to be modified by decreasing the number of low-et^Me clusters. This will 
also change the intrinsic distribution slightly. In the case of a purely oblate distribution, the 
opposite holds. An oblate cluster has a larger probability of being detected if its short axis 
is in the plane of the sky. This produces a bias towards more elongated clusters and one has 
to modify the distribution over projected ellipticities by increasing the number of low-et^Me 
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Figure 16. Relation between the ellipticity of a cluster generated in an Q = 1.0 universe and a cluster generated in an Q = 0.2 
universe, using the same random number seed to generate the clusters. For all nine clusters the ellipticities are calculated for 
the xy-, XZ-, and yz-projection, and using the tensor method of RHK. 

clusters. For the intermediate triaxial clusters (T = 0.5), the bias is much less, and the etme 
distribution has to modified only slightly. 



7 N-BODY SIMULATIONS 

As mentioned in the Introduction, the shapes of clusters will most certainly be infiuenced by 
the density parameter 0. For low 0, structures will collapse at a higher redshift than in an 
= 1.0 universe (Maoz 1990) and thus have had more time to virialize and wipe out their 
asymmetries and substructures. Furthermore, the timescale in a low-0 universe is longer 
than in an Finstein-de Sitter universe (see e.g. Padmanabhan 1993, eq. [2.79]). What may 
be also important is that in low-0 case, there is only little secondary infall at late times. 
This means that perturbations in a low universe cease to grow. On the contrary, = 1.0 
clusters keep accreting matter and this may change their structure considerably. 

In order to try and put constraints on 0, we have used N-body simulations of clusters in a 
CDM universe (van Kampen 1994). These simulations use a volume of radius 25 — 32h~^ Mpc, 
which is large enough to contain the turn-around radius and to model tidal fields. At a certain 
epoch galaxies are identified using a local density percolation technique combined with a 
virial equilibrium condition. For each galaxy found its constituent particles are replaced by a 
single softened particle, where the binding energy of the original particles is transferred into 
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Figure 17. Distribution over intrinsic ellipticities e^J^,t,l (solid line) and e^J^,t,2 (dotted line) for the simulated clusters with 
Q = 1.0. The ellipticities are calculated with the tensor method of RHK. 

Table 3. Mean ellipticities of the three orthogonal projections of the nine clusters, generated both using Q = 1.0 and Q = 0.2, 
and calculated with the three methods. The numbers in parentheses give the dispersions in e^. 

= 1.0 = 0.2 

(ei) 0.27 (0.01) 0.22 (0.01) 

(£2) 0.40 (0.01) 0.33 (0.02) 

(es) 0.26 (0.04) 0.29 (0.02) 

internal energy of the single particle galaxy. That the galaxy identification algorithm can 
be important was shown by Summers (1993). The projected density of the clusters on the 
'sky' satisfies that of a modified Hubble-profile (eq. [10]), and also the number of galaxies is 
comparable to the original data. For a complete description of the models and method see 
van Kampen (1994). 

For = 1.0, 75 clusters are generated according to the probability distribution of the 
parameters x, e and p, as given by Bardeen et al. (1986). These parameters determine 
the distribution of heights (;/), curvature (x) and shapes (e and p) of the peaks in the initial 
density field, out of which structures are assumed to form. A biasing parameter b = 2.0 
was adopted for the = 1.0 catalogue. Using a 2-D identification procedure, 4 out of the 
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Figure 18. Distribution over the triaxiality parameter T for the N-body simulated clusters with Q = 1.0. The ellipticities are 
taken from Figure 16. 

75 clusters are found to produce Abell-like clusters of richness class 0. All others model 
clusters resemble Abell clusters of richness R > 1. For this choice of the bias parameter the 
catalogue of simulated clusters is thus expected to be a fair representation of real clusters 
in the universe. Then each cluster is viewed from a random direction. All galaxies in the 
simulation cylinder that are projected onto the cluster area are used in order to include some 
background and foreground. For the two methods of RHK, we only use the galaxies inside a 
projected distance of lh~^ Mpc of the cluster center. This center is determined in the same 
way as was done for the real data. 

The distribution over uncorrected (=apparent) ellipticities is plotted in Figure 15. As 
before, the solid lines indicate the moments method of RHK, the dotted lines denote the 
tensor method of RHK, and the dashed lines are the results for the tensor method of PBF. 
This figure can be compared with Figure 1. It is obvious that the ellipticities of the simulated 
clusters are too large to be consistent with the observations. To examine the infiuence of 
on the shapes of the clusters, we have selected nine clusters (both massive and less massive) 
from the = 1.0 catalogue. These, then, are again simulated in an unbiased = 0.2 universe 
using the same random number seeds. The nine = 0.2 clusters are a representative sample 
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of a catalogue being constructed at present, analogous to the catalogue of = 1.0 clusters. 
The properties of this whole catalogue will be presented elsewhere. For each of these eighteen 
clusters (i.e. nine clusters both in = 1.0 and = 0.2) we have calculated €2 for the three 
orthogonal projections and using the tensor method of RHK (method 2). The mean values 
of ^-apparent for all three methods are shown in Table 3. The numbers in parentheses give the 
dispersions. The relations between €2 in the = 1.0 and the = 0.2 simulations are shown 
in Figure 16. From this figure it is clear that clusters in an = 0.2 universe tend to be more 
spherical than those in an = 1.0 universe. Fxamining all clusters individually, it turns 
out that almost all of the data having €apparent{^ = 0.2) > €apparent{^ = 1-0) correspond to 
clusters that either have just merged in the = 1.0 simulation and therefore appear more 
spherical in this simulation, or that are double clusters or appear to be due to the projection 
of two groups onto each other. All 'isolated' clusters turn out to lie below the dotted line. 
When the whole = 0.2 catalogue will be finished, we will have a handle on the fraction of 
clusters that are merging at present or that are/appear double. With this we will be able to 
put more quantitative constraints on the effect visible in Figure 16 and therefore on 0. 

Bardeen et al. (1986) show that the shapes of the density perturbations in the initial 
density field, which we expect to correspond with what we now call clusters, are somewhat 
more elongated for lower models. That their shapes are now less elongated points at a 
quite different evolution of clusters for different 0. For low 0, clusters collapse at a larger 
redshift, time elapses faster with redshift, and clusters do not accrete much material from 
outside during the last epochs. The latter means that = 0.2 clusters nearly have not grown 
during the last epochs and therefore their virialization process is not disturbed by clumps 
of infalling matter. 

It should be noted that these simulations show the same features as the original RHK 
data in that a cluster looks more elongated when considering only the brighter (i.e. more 
massive) galaxies. As mentioned previously, this is just an effect of the methods, and is 
not likely to be real. Also, the PBF method rejects a fraction of the clusters. Furthermore, 
clusters containing more galaxies are also intrinsically more spherical, though the effect is 
less pronounced than in the real data. 

To see how reliable it is to assume all clusters to have the same triaxiality parameter 
T (as assumed in Section 6), we determine the two (apparent) intrinsic ellipticities €int,i 
and €int,2 for the complete = 1.0 catalogue. The ellipticities are calculated using the 
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three-dimensional analog of the tensor method of RHK and the calculations are done within 
a (three-dimensional) radius of lh~^ Mpc, to resemble as good as possible the projected 
distribution. The results are shown in Figure 17. The solid line denotes the distribution over 
€int,i^ and the dotted line indicates €int,2- The distribution over T is shown in Figure 18. It 
is clear that more than 50 % of the clusters have a triaxiality parameter T larger than 0.8, 
which means that they are nearly prolate. Almost no clusters have T < 0.4 , which means 
that there are almost no oblate clusters in this sample. This result, although determined 
for the = 1.0 clusters and using the apparent ellipticities, is fully consistent with our 
analysis of the real clusters, where we found that the distribution over true ellipticities is 
better matched by a purely prolate distribution than a purely oblate one. It is not expected 
that this distribution will change drastically if we use the true intrinsic ellipticities instead 
of the apparent ones, because both €int,i and €int,2 will then change, and the effect on T may 
be small. 

8 DISCUSSION AND CONCLUSIONS 

We have reananlyzed a data set of 99 Abell clusters data, studied previously by Rhee, van 
Haarlem and Katgert (1989). The clusters are of Abell richness i? > 1, redshift z < 0.1, 
lO'' < a < 18^^, S < —25° and b > 30°. This data set is complete up to 2; = 0.08. To determine 
the shapes of the clusters, the same methods that these authors applied are used. Also, a 
method that has been used by Plionis, Barrow & Frenk (1991) in their analysis of the Lick 
catalogue, is applied. However, these three methods are found to have a number of errors and 
biases associated with them. By simulating clusters with a Monte-Carlo technique (using a 
Hubble-profile with a core radius of 0.25h~^ Mpc), these errors and biases are investigated. 
In fact, the adopted core radius is found not have a large influence in the range 0.1 — 0.5h~^ 
Mpc. Also the exact number of background galaxies adopted is of little influence. 

It is found that for small (projected) ellipticity etme (i-e., large axis ratio), the methods 
overestimate the ellipticity. This is due to the asymmetrical properties of the problem and 
is analogue to the Gott & Thuan effect (Thuan & Gott 1977). For large etme^ the methods 
underestimate its value due to the background and the spherical aperture (in the case of the 
Rhee et al. methods). 

The observed distributions over cluster shapes are corrected for the errors and biases, 
using the iterative Lucy-algorithm. The corrected distribution shows a peak at etme ~ 0.4 
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and extends to ~ 0.8. The number of nearly spherical clusters (with etme ^ 0.2 ) is limited 
to about 20%. These results are consistent with some of the previous studies, but have much 
better statistics due to both the larger number of clusters in our sample and the larger 
number of galaxies per cluster. Furthermore, the result is internally consistent. That is, we 
are able to reproduce both the observed distributions over tapparent and the observed relation 
between tapparent and the number of galaxies in a cluster. This is only possible if one assumes 
an anti-correlation between the number of galaxies in a cluster and its true ellipticity. Using 
redshift-independent richnesses of Vink & Katgert (1994), defined as the number of galaxies 
in the absolute magnitude interval [-22.5,-21.5], we show that the number of galaxies in a 
clusters is correlated with the richness of that cluster and is not just an effect of the cluster 
distance. This means that richer clusters are intrinsically more spherical than poor ones. It 
may point to the fact that rich clusters have turned around from the Hubble-fiow earlier 
and have had more time to virialize. Also, the virialization process will be faster due to the 
larger number of galaxies. To achieve the internal consistency, it is not necessary to assume 
that clusters are more elongated if only the brighter galaxies are included (as was suggested 
by Binney 1977). This turns out to be only an artifact due to the smaller number of galaxies 
used to determine tapparent- 

We try to say something about the intrinsic shapes of clusters by assuming that they 
all have the same triaxiality parameter T, as defined by Franx, lUingworth & de Zeeuw 
(1991). If one assumes all clusters to be oblate (i.e., have T = 0.0), most of them will have 
intrinsic ellipticities of ~ 0.55. In this case almost no intrinsically spherical systems are 
found and a rather large percentage of highly fiattened systems is predicted. In the case 
of a purely prolate distribution (T = 1.0), the average intrinsic ellipticity is also ~ 0.55. 
However, slightly more clusters with small elongations are present than in the oblate case, 
and less with high elongations. When only intermediate triaxial clusters (having T = 0.5) 
are considered, the profile is intermediate between those for oblate and prolate shapes. If 
we try to reconstruct the profile over etme from these constant-T distributions, the T = 
0.5 and T = 1.0 distributions perform best. The purely oblate population produces too 
many spherical clusters. This result is in agreement with the analysis of Salvador-Sole & 
Solanes (1993), who find the projected elongations of clusters to be consistent with a prolate 
distribution with intrinsic axial ratios Gaussian-distributed with mean ~ 0.5 and standard 
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deviation ~ 0.15. These elongations are mainly produced by the tidal interactions of massive 
enough nearby clusters. 

Because the density parameter will most certainly have an influence on the shapes 
of clusters, we use N-body simulations, which include a recipe for galaxy formation and 
merging. A CDM-spectrum and a density parameter = 1.0 are used. The biasing parameter 
found for these simulations is b = 2.0. The simulations form a complete catalogue that is 
expected to give a fair representation of all real clusters. The projected galaxy distributions 
follows roughly the projected Hubble-proflle, eq. (10). 'Observing' these clusters in the same 
way as the real clusters produces too large elongations as compared to the observations. For 
= 0.2 a similar catalogue is being prepared, though not ready yet. We have picked out 
nine of the = 0.2 clusters (both massive and less massive ones) that are simulated in the 
same way as their = 1.0 counterparts, using the same random number seeds. A biasing 
parameter 6 = 1.0 is adopted here. The = 0.2 clusters are, on average, less elongated than 
their = 1.0 counterparts. The cases where the = 0.2 cluster is more elongated turn out to 
be due to merging, projection of two clusters onto each other, or because the = 0.2 cluster 
actually is a double cluster. The preliminary results of the N-body simulations suggest that 
= 0.2 is more consistent with the data than for = 1.0. A more quantitative analysis will 
be possible when the whole = 0.2 catalogue will be flnished. The results of the simulations 
are consistent with the results of Evrard (1993) and Evrard et al. (1993), who use simulations 
of the X-rays tracing the cluster potential. They also found low-0 clusters to be more nearly 
spherical and also more centrally condensed than = 1.0 clusters. However, Evrard (1993) 
concludes that their = 1.0 simulations give a better description of the observations than 
do the = 0.2 simulations, although they have only 4 clusters to compare with. A follow-up 
paper with a sample of 50-60 EINSTEIN-clusteis is in progress (Mohr et al. 1994). 

As the = 0.2 catalogue is not complete yet, we use the = 1.0 catalogue to calculate 
the intrinsic shapes of the model clusters. They turn out to be mainly nearly prolate, with 
more than 50% having a triaxiality parameter larger than 0.8. Almost no (nearly) oblate 
clusters are detected. This is consistent with what is expected from the distribution of true 
projected cluster shapes. 
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